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ABSTRACT 

For a wide variety of initial and boundary conditions, adiabatic one dimensional flows of an ideal gas 
approach self-similar behavior when the characteristic length scale over which the flow takes place, 
i?, diverges or tends to zero. It is commonly assumed that self-similarity is approached since in the 
R oo(0) limit the flow becomes independent of any characteristic length or time scales. In this 
case, the flow fields f{r,t) must be of the form f{r,t) = t"f F{r / R) with R cx (±t)". We show 
that requiring the asymptotic flow to be independent only of characteristic length scales implies a 
more general form of self-similar solutions, f{r,t) = R^fF{r/R) with R cx R^ , which includes the 
exponential {S = 1) solutions, R oc e*/'^. We demonstrate that the latter, less restrictive, requirement 
is the physically relevant one by showing that the asymptotic behavior of accelerating blast-waves, 
driven by the release of energy at the center of a cold gas sphere of initial density p oc r~", changes 
its character at large uj: The flow is described by < (5 < 1, i? oc t^/'-^~^\ solutions for uj < uJc, by 
(5 > 1 solutions with R oc (— diverging at finite time {t = 0) for uj > uJc, and by exponential 
solutions for lu = uic {uJc depends on the adiabatic index of the gas, cjc ^ 8 for 4/3 < 7 < 5/3). 
The properties of the new solutions obtained here for tu > lUc are analyzed, and self-similar solutions 
describing the t > behavior for uj > ujc are also derived. 
Subject headings: hydrodynamics — shock waves — supernovae: general 



1. INTRODUCTION 

Self-similar solutions to the hydrodynamic equations 
describing adiabatic one dimensional flows of an ideal 
gas are of interest for several reasons. The non-linear 
partial differential hydrodynamic equations are reduced 
for self-similar flows to ordinary differential equations, 
which greatly simplifies the mathematical problem of 
solving the equations and in certain cases allows one 
to find analytic solutions. Moreover, self-similar solu- 
tions often describe the limiting behavior approached 
asymptotically by flows which take place over a char- 
acteristic scale, R, which diverges or tends to zero (see 
lSedov|[l959t iZel'Dovich fc RaizedfTOGTt lBarenblattl[l996l 
for reviews). Some examples of such asymptotic solu- 
tions which are widely used in astrophysical contex ts are 
the Sedov-voii Neumann- Taylor solutions (Sedovl IT946I : 
I Von Neuma^ll947l: lTavlodfl950l) describing expanding 
decelerating spherical blast waves, for which R ^ 00, 
and t he Gandel'Man-Frank-Kamenetskii-Sakurai solu- 
tions (iGandel'Man fc Frank-Kamenetskiilll956t ISakurail 
119601 ) describing the emergence of a shock wave from the 
surface of a star, for which i? — )■ 0. Both types of solu- 
tions are relevant, e.g., to supernova explosions and in 
particular to the recently detected shock breakouts (e.g. 
iMatzner fc McKed [l999l : IWaxman et al.ll2007t ). An ex- 
tensive discussion of spherical self-sir nilar blast-waves in 
an astrophysical context is given by lOstriker fc McKei 
( fr98S) . 

It is commonly assumed that self-similarity is ap- 
proached since in the R oo(0) limit the flow becomes 
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independent of any characteristic length or time scales. 
Using dimensional arguments it is possible to show that 
if the flow is determined by a set of constants, using 
which it is impossible to construct a constant with the 
dimensions of length or time, then the flow fields f{r,t) 
must be of the for m f(r,t) = t°'-f F(r/R) with R oc (±t)" 
(see chapter XII of IZel'Dovich fc Raizeijll967( ). We show 
in § [2] that requiring the asymptotic flow to be inde- 
pendent only of characteristic length scales is sufficient 
for showing, based on dimensional arguments, that the 
flow must be self-similar. The less restrictive require- 
ment allows a more general form of self-similar solutions, 
f{r,t) = R^fF{r/R) with R oc R^, which includes the 
exponential (5 = 1) solutions, R oc e*/"^. The existence 
of exponential self-si milar solutions has been noted by 
several authors (e.g. IStanvukovichI |1960() . However, it 
was generally assumed that asymptotic solutions, which 
are of interest, are of a power law form. 

In § [3] we show that the asymptotic self-similar so- 
lutions describing the propagation of accelerating blast 
waves, propagating in a cold gas sphere of initial den- 
sity p oc r~" with a; > 3, are of the more general 
form, R (X R^ , with exponential solutions obtained at 
UJ = ujc{"f), where 7 is the adiabatic index of the gas. 
The new solutions obtained here for uj > ujc extend the 
family of second type solutions describing the asymp- 
to tic flow of acc elerating blast waves, which was derived 
bv lWaxman fc Shvarts, (1993. ) and was limited to uj < uJc, 
to UJ > UJc- The properties of the new solutions are an- 
alyzed in § 13.21 and self similar solutions describing the 
UJ > UJc flow at times later than the finite divergence time 
are derived in § 13.31 Our results are summarized in § |4l 

It should be noted here that the self-s imilar solu- 
tions derived by IWaxman fc ShvartsI ()1993| ) exist only 
for UJ > ujg{'-f) > 3, where ujg = 3.26 for 7 = 5/3 and 
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approaches 3 for 7 — >■ 1, while the Sedov-von Neuraann- 
Taylor solutions provide the correct asymptotic solutions 
only for a; < 3. The asymptotic behavior within the (nar- 
row) range of 3 < w < ^^(7) is not described by either 
of the two types of solutions. The nat ure of the asym^ 
totic flow in this regime is discussed in iGruzinovl (|200 



totic now m tins regime is di 
iKushnir fc WaxmanI ()2010[ ). 



2. THE GENERAL FORM OF ASYMPTOTIC 
SELF-SIMILAR SOLUTIONS 

The equations describing adiabatic one-dimensiona l 
flow of an ideal gas are (e.g. IZel'Dovich fc RaizedfTOGl 

{dt+udr)u + p-^dr{-f~^pc'^)^0, 

idt+udr){c^p'-')^0, (1) 

where u, c, and p are the fluid velocity, sound speed and 
density respectively (the pressure is given hy p = pt? 
and V = 1,2,3 for planar, cylindrical and spherical sym- 
metry respectively. 

Consider a solution /(r, i; {cj}, 7), where / stands for 
p, u or c, and {ci} are the set of constants determin- 
ing the initial and boundary conditions. Let us as- 
sume that the flow takes place over a characteristic 
scale which diverges or tends to zero monotoni- 

cally with time (if R{t) is not monotonic then the flow 
is characterized by some finite length scales correspond- 
ing to the extrema of R). In this case we may replace 
the variables {r,t\ with {^(r, i) = r/R,R}, and de- 
scribe the flow by /(^, R; {a}, 7), where f{r, t; {a}, 7) = 
f[^{r,t),R{t); {c,;},7], and by R{R; {ci},-f), which deter- 
mines the relation between R and t. 

Let us assume next that in the limit R — >■ oo(0) the 
flow becomes independent of any characteristic length 
scales. That is, that in the limit R — >■ oo(0) the flow 
is determined by only a subset of {ci}, denoted {bi}, 
from which a constant with the dimensions of length 
cannot be constructed. The dimensions of each param- 
eter bk may be considered as a three dimensional vec- 
tor Vk, [bk] = M"lL^=^T^^ where M, L, and T are the 
units of mass, length, and time, respectively. Since a 
constant with the dimensions of length may not be con- 
structed using there are at most two constants 
with independent dimension vectors. As we shall see 
below, there must exit two constants with independent 
dimensions. Denoting these by {01,02}, and express- 
ing the other constants as bi = 7^0^" 02'^ with di- 
mensionless 7^, the asymptotic solution may be written 

as /as.(?, ^;ai,a2,{7j;7), Ras.{R;ai,a2,{'ji},j). Since 
the dimensions of {/?, 01,02} are independent (otherwise 
it would be possible to construct a constant with dimen- 
sions of length using {oi, 02}), there is a single combina- 
tion (product of powers) of {i?, 01,02} that has a given 
dimension vector. This implies that /as. must be of the 
form /as. = o"^'o2^^i?"^i^(^; {7i},7) and i?as. must be 

of the form 7?as. = aN2'i?*i^i?({7j,7)- 

Thus, the assumption that the asymptotic solution is 
independent of any characteristic length scales implies 
that it must be a self-similar solution of the form 

u^mUiO, c = mC{0, p = BR'G{0 (2) 
R = AR^. (3) 



with 



The spatial part of u{c) was chosen as ^t7(^) (^C(^)) for 
convenience (see below). Note, that the dimensions of 
A and B are independent. These solutions include the 
R oc (±t)" solutions, obtained for 5 = (a — l)/a ^ 1, and 
the exponential R oc e*/"^ solutions obtained for ^ = 1. 

Substituting the ansatz of eqs. ^ and ^ in the hy- 
drodynamic equations (eq. ([T])), the partial differential 
equations are replaced with a single ordinary differential 
equation, 

df/ _ Ai({/,C) 



dC A2(C/,C)' 



and one quadrature 
dXxvS, A{U,C) 



dU Ai(C/,C) 
G is given implicitly by 

(^C)-2(''+^)|l _ C/|Ac.(7-l)(^+e)+A^<>A ^ CoilSt., (6) 

with 

A = -(7- l)e-f 25. (7) 
The functions A, Ai, and A2 are 

A=c2 - (1 -c/)^ 

£ + 2(5 



Ai ^U{l-U){l-U -S) - C l^i^U + 
A2^C{{1-U){1-U -S) 

^l-lu[ii.-l){l-U) + S]-C' 
25-(7-l)e C2 



27 



1 -U' 



(8) 



3. ACCELERATING BLAST WAVE SOLUTIONS 



Consider the blast wave produced by the deposition of 
energy E within a region of characteristic size d at the 
center of an initially cold (p = at r > d) gas sphere 
with initial density po = Kr~'^ (at r > d). As the shock 
radius R diverges, it is reasonable to assume that the flow 
becomes independent of any length scale determined by 
the initial conditions, i.e. that R{t) is the only relevant 
length scale. In this case, the flow should approach a 
self-similar solution of the form given by eqs. ^ and ([31). 
Since the density just behind the shock wave is, for strong 
shocks, a constant factor, (7 -|- l)/(7 — 1), times the den- 
sity just ahead of the shock, we must have e = — w, and 
we may choose B = K. With this normalization, the 
Rankine-Hugoniot relations at the shock front determine 
the bou ndary conditions for the se lf-similar solutions to 
be fe.g. IZel'Dovich fc Raizeil[l967l ) 



f/(l) 



7 + 1' 



C(l) 



v/27(7 - 1) 
7+1 '' 



G(l) 



7 + 1 



(9) 

The only parameter of the self-similar solution that re- 
mains to be determined is 6. The methods described 
below for determining S, and for analyzing the result- 
ing solutions^roperties, are similar to those described 
in (jWaxman fc Shvarts 1993). The discussion in the fol- 
lowing subsections is therefore very con cise. The reader 
may refer to iWaxman fc Shvarti ()1993D for more elabo- 
rate explanations. 
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3.1. U! < UJc 

In the Scdov-von Neumann- Taylor analysis it is as- 
sumed that the second dimensional constant, in ad- 
dition to K, that determines the self-similar solution 
is E. In this case dimensional considerations imply 
R oc (Et^/K )^/^^-'^\ i.e. 6 = (to - 3)/ 2. As explained 
in detail in (IWaxman fc ShvartsI 1 1 9931 ). the Sedov-von 
Neumann- Taylor solutions are the correct asymptotic so- 
lutions only for w < 3, for which S < and the blast wave 
decelerates with time. For larger values of ut the mass 
and energy contained within the self-similar solution are 
infinite, reflecting the fact that the initial gas mass at 
r > d diverges for d 0. It was therefore suggested by 
IWaxman fc ShvartsI (|1993D that for a; > 3 the asymptotic 
solution is given by a self-similar solution only over part 
of the {£,,R} plane, bounded by ^ = 1 and £,c{R) < 1, 
and by a different solution at < ^ < £,c{R)- 

Since ^c(-R) must be a contact or a weak discontinuity, 
£,c{R) must be a characteristic of the self-similar solution. 
For the self-similar flow, the characteristic lines 

ri ^''0 ^ dr± 

Co : — = M, C±: —— = M ± c, (10) 
dt dt 

are given by 

The flow just behind the shock is always subsonic: the 
shock-front point {U = U{1),C = C(l)} lies above the 
"sonic-line" U + C = 1, which implies that C+ char- 
acteristics emerging from points just behind the shock 
always overtake it. C+ characteristics that do not over- 
take the shock exit only if the self-similar solution crosses 
the U + C = 1 line in the {U, C} plane into the region 
where U + C < 1. Requiring to coincide with a C+ 
characteristic, that does not overtake the shock, implies 
therefore that the solution must cross the U + C = 1 line. 
Since A = along the sonic line, eqs. ([5]) imply that a 
physical solution must cross the sonic line at a singular 
point Ai = A2 = (otherwise U{£,) or C(^) are not 
single valued). This requirement determines the correct 
value of S for the w > 3 asymptotic solutions. 

The self-similar solutions o btained in this way for ui > 
3 were analyzed in detail in (jWaxman fc ShvartsI 1 19931 ). 
They describe accelerating blast waves, with 5 > 0, and 
approach the singular point {U = 1 — 5, C = 0} as 
^ — >■ 0. Analyzing the behavior of the solutions near 
this singular point, it was shown that although the mass 
and energy contained in the self-similar solution are in- 
finite, the mass and energy contained within the region 
£,c{R) < ^ < 1, where £,c{R) is a C+ characteristic which 
satisfies ^c(-R) ^> as i? — >■ 00, approach finite values as 
i? — !• 00. Moreover, it was shown that £,c{R)R oc t as 
i? — 00, implying that the asymptotic flow within the 
region < r < £^c{R)R is described by the self-similar so- 
lution of e xpansion into vacuum. Fin ally, it was demon- 
strated in (jWaxman fc Shvartslll993l) by numerical sim- 
ulations that the asymptotic behavior described above is 
indeed approached for R/d^ 1. 

3.2. uj > UJc 



1.5 1 ! 1 1 1 1 r 
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Fig. 1. — S{ui) for 7 = 5/3 (red dashed line) and 7 = 4/3 
(solid black line) determined by the requirement that the C(U) 
curve crosses the sonic li ne at a singular point. As noted in 
llWaxman fc ShvartsI IT993) . there are no solutions satisfying this 
requirement within some narrow range 3 < a; < u]g{f) {uig = 3.26 
for 7 = 5/3 and approaches 3 for 7 — >■ 1). Th e nature of the 
asym ptotic flow in this regime is discussed in IIKushnir fc WaxmanI 
[20101) . 



In the analysis of ("Wax man fc ShvartsI Il993l ) it was 
assumed that the asymptotic solutions are of the form 
R cc t". It was found that a diverges as uj approaches 
a finite value, 1^-1^(7) (wc ~ 8 for 4/3 < 7 < 5/3 a nd di- 
verges for 7 — > 1, see fig. 5 of IWaxman fc Shvart!j ri993'). 
Allowing solutions of the more general form, i? oc i?* 
as suggested in § [21 the divergence of a suggests that 
the asymptotic solution for lo = uJcil) is an exponential 
solution, R oc e*/'^ (recall that 6 ~ (a — l)/a). Search- 
ing for solutions of the form R (x R^ , that cross the 
sonic line at a singular point, we find that such solutions 
do exist for uj > lUc- i5(w) is shown for 7 = 5/3 and 
7 = 4/3 in figure [TJ We find 5 = 1 for cj = ujc{'^), im- 
plying exponential expansion, and S > 1 for uj > Wc(7)i 
implying that the shock radius diverges in finite time, 
R oc {—ty/'--^~^\ Fig. 12] presents a comparison of the 
C{U) curves obtained by numerically solving the hydro- 
dynamic equations, eq. ([T]), with those obtained by solv- 
ing the self-similar eq. (U). The numerical results suggest 
that the self-similar solutions provide a correct descrip- 
tion of the asymptotic flow. 

For (5 > 1, the self-similar solution approaches the sin- 
gular point {U — 0,C = 0} as^ 0. Solving the 
self-similar equations in the vicinity of this point we find 
that in the limit of ^ the solutions are given by 

C/ = --^, C oc C/(3^-i)/2, Goce""(lnO""^ (12) 

for (5 = 1, and by 

u = Uie.^-^, C = Cie'^-\ G = Gir", (13) 

where J7i, Ci, Gi are constants, for (5 > 1. Using eqs. 
and (|13p we find that the G+ characteristics originating 
from points below the sonic line are given in the limit 

i? —5- 00 by 
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-uj = 7.686, 5 = 1.0, Self-Similar 
-oj = 7.686, Numeric 
-to = 13.0, 5 = 2.1, Self-Similar 
-uj = 13.0, Numeric 
-Sonic Line 




Fig. 2. — A comparison of the C(U) curves obtained from numer- 
ical solutions of the hydrodynamic equations, eq. JTJ, with those 
obtained by solving the self-similar eq. Q. Results are shown for 
7 = 5/3 and two values of oj, oj = 7.686 for which the self-similar 
analysis yields <5 = 1.0 (lower curve), and oj = 13.0, for which the 
self-similar analysis gives 5 = 2.10. Numerical solutions are shown 
for R/d = 500 (where R is the shock radius and d the radius of the 
sphere within which the energy is initially deposited). 



The energy contained in a region of the self-similar so- 
lution, bounded by ^ = 1 and a C+ characteristic £,+ {R), 
is given by 

\2 ij3-w+2i5 



E+{R) =AtiKA^R 

X / d^^^G 



7(7 - 1) 



,(15) 



In order for the energy not to diverge as i? — > cxd and 
£,+ {R) — > 0, we must have 



S < 



(16) 



This is satisfied for the solutions presented in fig. [TJ For 
S values satisfying eq. (jl6p . the integral over ^ diverges 
as S,+ {R) 0, while i?3-w-f25 vanishes as i? ^ 00. Us- 
ing eqs. (|12|) -(fT4 |) . it is straight forward to verify that 
E+{R) Const, as i? 00. Thus, the energy con- 
tained in the region of the flow which is described by 
the self-similar solution, £,c{R) < ^ < 1, approaches a 
finite value as i? 00. A similar calculation shows that 
the mass contained in this region also approaches a finite 
value. 

3.3. t > 

For S > 1, the shock radius diverges at a finite time, t = 
0. This implies that the spatial distribution of the fiow 
fields at the divergence time is described by the ^ — > 
behavior of the solution, which is given by eq. (fT3l) . That 
is, at i = the spatial distribution of the flow fields is 

uocr'', c (X u, p oc r^" (17) 

(Note, that this form of the flow flelds applies for r > vq 
where rg is the finite radius approached by ^cR as R 
diverges, see eq. (fT4|l ). Since the distribution of flow 
fields at i = has no characteristic scale, the flow at 
t > will be described by a self-similar solution. 



The t > self-similar solution may be constructed as 
follows. Denoting the scaling radius of this solution by 
R, we choose 



R = -AR" 



(18) 



with the same values of A and S obtained for the i < 
solution. This implies R oc t^^'^^^^), i.e. R is decreasing 
with t and diverges a,s t ^ 0+. Next, we choose a C{U) 
curve that approaches in the limit ^ ^ 0, where ^ = r/ R, 
the same singular point approached in this limit by the 
t < solution, [/ = C = 0: 



■<5-l 



c = Cie 



<5-l 



G = GiC 



(19) 



with the same values of Ui,Ci, Gi obtained for the t < 
solution (Note, that for (5 > 1 any straight line C oc U 
is a solution of eq. (j4|) near U — C — 0). It is straight 
forward to verify that the t > self-similar flow thus 
deflned approaches the form given by eq. pT)) as < — > 0+ . 

In the limit t ^ 00, R ^ and the spatial distribu- 
tion of the flow flelds is described by the f ^ 00 be- 
havior of the solution. In the limit ^ — )■ 00, the C{U) 
curve of the t > solution approaches the singular point 
{U = 1 — 6,C = 0}. The behavior of the self-similar 
solutions in the vicinity o f this point was analyzed in 
(jWaxman fc ShvartdfToM . 

4. SUMMARY 

We have discussed the asymptotic behavior of adia- 
batic one dimensional flows of an ideal gas, which take 
place over a characteristic scale R that diverges or tends 
to zero. We have shown in §[2] that requiring the asymp- 
totic flow to be independent of characteristic length 
scales implies, based on dimensional arguments, that the 
flow must be self-similar, with flow fields given by eq. ^ 
and R satisfying eq. R oc R^ . The ordinary differ- 
ential equations determining the self-similar profiles are 
given by eqs. g])-®. 

In § [3] we have shown that the asymptotic self-similar 
solutions describing the propagation of accelerating blast 
waves, propagating in a cold gas sphere of initial density 
p oc r~" with w > 3, are of the general form R oc R^ , with 
exponential solutions obtained at w = u)ci"f)- <5(^^,7) is 
shown in fig. ^ for 7 = 4/3,5/3. Fig. [2] demonstrates 
that numerical solutions of the hydrodynamic equations, 
eq. (HI) , indeed approach these self-similar solutions as R 
diverges. The properties of the self-similar solutions ob- 
tained for UJ > u>c{j) are analyzed in § 13.21 The C{U) 
curves of these solutions approach the singular point 
{U = 0, C = 0} as ^ — 0. Analyzing the behavior of the 
solutions near this singular point we have shown that the 
energy and mass contained in region of the R} plane 
bounded by ^ = 1 and ^ — ^+{R), where ^+{R) is a C+ 
characteristic that approaches ^ = as i? diverges, both 
tend to finite constants as R diverges. 

For 5 > 1, the shock radius diverges at a finite time, 
t = 0. This implies that the spatial distribution of the 
flow fields at the divergence time is described by the ^ — ?> 
behavior of the solution, which is given by eq. ([T^ . 
This implies that at i = the spatial distribution of the 
flow flelds is given by eq. (I17|) . In § 13.31 we have shown 
that the t > flow is described in this case by a self- 
similar solution with the same value of S as that of the 
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t < solution, see eq. 
by eq. ([T9| . 



181 , and spatial profiles determine 
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